Electronic Mechanism of Superconductivity in the Cuprates, Ceo? and Polyacenes 
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On the basis of an analysis of theoretical and numerical studies of model systems, and of experiments 
on superconductivity in doped Ceo, polyacenes, and the cuprate high temperature superconductors, 
we propose that a purely electronic mechanism of superconductivity requires structures at an in- 
termediate or mesoscale. Specifically, we address the crucial question of how high temperature 
superconducting pairing on the mesoscale can arise from purely repulsive electronic interactions. 

PACS: 

I. INTRODUCTION 

Induced attractions between electrons, such as are me- 
diated by phonons in conventional superconductors, are 
typically weak while the direct repulsive interactions be- 
tween electrons are strong. It has therefore long been 
felt that the key to high temperature superconductivity 
is a purely electronic mechanism that directly exploits 
the repulsive interaction to produce pairing. 

After the discovery of high temperature superconduc- 
tivity in the cuprates by Bednorz and Miiller an idea 
for how this could be achieved was proposed 1|. It was 
suggested that an antiferromagnet disordered by quan- 
tum fluctuations and characterized by resonating singlet 
pairs of spins could be viewed as a state in which super- 
conducting pairing occurs, even though the system is an 
insulator with one electron per site and zero superfluid 
density! In this case, when doped, so that charge mo- 
tion is allowed, the system will automatically become a 
superconductor. A related idea is that the singlet pair 
formation in the insulating state produces a gap in the 
spin excitation spectra HU, which evolves into the su- 
perconducting pairing gap upon doping. 

A fundamental flaw with this approach soon became 
clear; in the relevant two-dimensional spin 1/2 antiferro- 
magnet on a square lattice, the ground state is magneti- 
cally ordered || , and indeed so strongly that the elemen- 
tary excitations about the ordered state, spin waves, give 
a quantitatively accurate description of many quantities. 
Needless to say, the magnetically ordered state bears no 
resemblance to a superconductor, even locally! 

Yet, the intuition that the high energy scale of high 
temperature superconductors must have its origin in 
the repulsive interactions remained. The discovery by 
Hebard et al. || of high temperature superconductiv- 
ity in alkali doped Cgo was another milestone. It illu- 
minated a path, in the clearest terms, that the missing 
element in the above revolutionary argument is the exis- 
tence of a mesoscale structure. The Ceo molecule itself 
plays the role of this structure permitting electronic pair- 
ing, which would disappear if the molecule were expanded 
to the macroscopic scale. We were thus led to propose 



a mechanism |7|-|9[ of superconductivity in which attrac- 
tion arises from repulsion. Later it was suggested that a 
similar pairing mechanism may apply in the cuprate high 
temperature superconductors 1 10| , where the relevant 
mesoscopic structures are self-organized stripes. 

Here, we present evidence, which we consider to be 
exceedingly strong, that this set of ideas underlies the 
occurrence of high temperature superconductivity in di- 
verse materials. It is shown that local singlet formation 
as a driving mechanism for "pairing" is correct, but only 
on a mesoscopic scale. The pairing tendencies are either 
much weaker or non-existent in extended systems. 

There are also a number of remarkable recent exper- 
iments on doped Cgo [ p^|JT^ ] and polyacene crystalline 
films |l4| which encourage the view that these are elec- 
tronically driven true high temperature superconductors 
in disguise. Although the authors of these papers ascribe 
the superconductivity of these materials to the electron- 
phonon interaction, we feel that the experimental facts, 
including the high value of T c itself (presently up to 52K), 
rather indicate that these materials are true "high tem- 
perature superconductors," in which correlation effects 
play a dominant role both in the normal state and in the 
mechanism of superconductivity. We shall briefly discuss 
this issue in the final section. 

It is important to stress the distinction between pos- 
sible and necessary consequences of the existence of 
mesoscale structures — surely not all forms of repul- 
sive interactions on any given cluster will lead to pair 
binding. Additionally, there are strong reasons [|l5|,[l6| to 
expect the same mesoscale physics to lead to competing 
orders. Indeed, we propose that the existence of a com- 
plex phase diagram with a variety of broken symmetry 
phases, including high temperature superconductivity, is 
a characteristic feature of a material with an electronic 
mechanism of pairing on the mesoscale. 



II. MESOSCALE PAIR-BINDING 

At atomic scales, the Hamiltonian is simple (although 
strongly interacting), and the dominant Coulomb inter- 
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actions are readily identified. At the macroscopic scale, 
the physics is again simple and is governed by the univer- 
sal properties of the attractive fixed points of the renor- 
malization group theory. However, as the high energy 
degrees of freedom are integrated out, beginning with 
the microscopic scale, every imaginable multiparticle in- 
teraction is generated, so that the effective Hamiltonian 
becomes complicated and the dominant physics can be 
obscure. Thus, it is necessary to first determine whether 
there is any robust and predictable behavior at all. We 
argue that on certain finite clusters quantum fluctuations 
of the spins can produce a disordered ground-state with 
a large gap in the spin excitation spectrum and a robust 
pair-binding of holes or electrons. 

In each case, we start with a "neutral" cluster with N 
sites and one electron per site. For all the systems stud- 
ied, the neutral cluster has a unique, spin singlet ground- 
state. We define E(Q, S) to be the lowest energy eigen- 
value of the cluster with total spin, S, and "charge," Q, 
or equivalently a total of N + Q electrons. Whenever the 
ground-state is a spin singlet, we can define the spin gap 
A s to be 



A S (Q) = E(Q,1)-E(Q,0). 



The pair-binding energy is defined as 

E P (Q) = 2E{Q + 1) - E(Q + 2) - E(Q) 



(1) 



(2) 



where E{Q) signifies the minimum over S of E{Q,S). 
Clearly, a positive pair-binding energy signifies an effec- 
tive attraction between electrons, in the sense that given 
2(N + Q + 1) electrons and two clusters, it is energetically 
preferable to place N + Q + 2 electrons on one cluster and 
N + Q on the other than to put N + Q + 1 electrons on 
each of the clusters. 

We illustrate the origin of these effects, the relation 
between A s and E pi and the importance of intermediate 
scales to the occurrence of such "attraction from repul- 
sion" by analyzing a variety of models. We defer to the 
following section any discussion of the relation between 
spin-gap formation and pair-binding with the occurrence 
of superconducting or other long-range order in extended 
systems. 



A. Hubbard rings 

As a first example consider Hubbard rings with N sites, 
where N is even. The Hamiltonian is 



(3) 



with the implied periodic boundary condition. The 
fermion operator c| CT creates an electron with spin a on 
site i, and rii a = c\ a Ci a is the density of electrons of spin 
a on site i. 



The energies, E(Q, S), can be calculated exactly from 
Bethe Ansatz |l7]] . We have extended existing results for 
the pair-binding energy, E p , by numerically evaluating 
the Bethe- Ansatz equations for a wide range of system 
sizes as large as 1024 sites, and by calculating the spin- 
gaps, A s . The exact particle-hole symmetry of the Hub- 
bard model on a bipartite lattice implies that electron 
doping, Q > 0, and hole doping, Q < 0, are equivalent. 
As can be seen in Fig. [I] pair-binding occurs for elec- 
trons added to the neutral molecule, E p (0) > 0, whenever 
N — An, but does not occur when N — An + 2, where n is 
a positive integer. The difference between these two cases 
can be readily understood from low order perturbation 
theory in U / 1 as discussed in Appendix |a|. 

The role of intermediate scales can be seen directly 
from Fig. |l|; the pair-binding energy vanishes for large N 
and is maximum at an intermediate value of N . More- 
over, in Fig. H we show E p for fixed TV as a function of U, 
from which it is clear that it is maximal for intermedi- 
ate interaction strengths as well. Indeed, we have proven 
a rigorous theorem |Q that for the standard Hubbard 
model on any finite lattice and with positive hopping 
matrix elements, E p < for U — > oo. (Trivially, E p < 
for U = 0.) 
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FIG. 1. Pair-binding energy, E v , of N = An and TV 
site Hubbard rings with t = 1 and (7 = 4. 
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The crucial question is what is the mechanism of pair- 
binding? We can gain some insight into this question 



in the large- N limit from bosonization |19|. We show 
that pair-binding is closely related to the phenomenon 
of spin-gap formation. It is well-known that at long 
wavelength, spin-charge separation occurs in ID Hub- 
bard rings. Therefore, all energies can be written as a 
sum of a spin contribution and a charge contribution. 
The result is that for N = An > 1, 
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FIG. 2. Pair-binding energy, E p (solid symbol), and 
spin-gap, A s (hollow symbol), of 12 site Hubbard ring as a 
function of U in units of t = 1. 
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Bx ln 1/2 (N) + B 2 



E n — A, + — — — — 



(4) 
(5) 



Here, v s and v c are the spin and charge velocities, respec- 
tively (in units in which the lattice constant is unity) and 
A c is the charge-gap in the N — > oo limit. The constants, 
Bj , are numbers of order unity, which we have not evalu- 
ated. The large- N behavior of these quantities computed 
from the Bethe Ansatz are shown in Fig. |^ and appears 
to be in good agreement with these expressions. 

Because of the relevance of Umklapp scattering, the 
charge degrees of freedom are described by a sine- Gordon 
field theory. Dilute doping results in a dilute gas of spin 
charge ±e solitons (holons or eons [pj) with a gapped 
spectrum. These solitons can, in turn, be mapped onto 
massive Dirac fermions. Consequently, the finite size 
level spacings are those of a non-relativistic particle in 
a box, which gives rise to the third term in Eq. How- 
ever, in making a spin excitation in the neutral molecule, 
no charge is added, so the spin-gap is independent of the 
charge dynamics. 

It has been previously observed |^0|,|l| that, because 
the spin-system is quantum critical, the spin correlation 
functions are conformally invariant, so that the spin-gap 
must scale with v s /N. The slow renormalization group 
flow associated with a marginally irrelevant operator at 
the fixed point invalidates this finite size scaling behavior, 
however. A consequence of this for the infinite system is 
that the spin-spin correlation function falls with the dis- 
tance, r, as (ln 1 ^ 2 r)/r p^] . We conjecture that a similar 
logarithmic correction occurs in the finite size scaling be- 
havior, hence the logarithmic term in Eq. (Q). As can 
be seen from the large N behavior of the Bethe Ansatz 
solution in Fig. |3|, just such a logarithmic term appears 
in the data. In terms of the elementary spinon excita- 



tions of the spin 1/2 Heisenberg model, A s is interpreted 
as the energy of a two-spinon state, while A s + B 3 v s /N 
is twice the single spinon energy; —B 3 v s /N thus has the 
interpretation of an effective spinon-spinon interaction. 

An important lesson of this analysis is that, for large 
N, E p w A s . This relation, originally proposed by us 
in the context of pair-binding on a Cqq molecule ||, is 
in this case a consequence of the quantum critical char- 
acter of the ID Hubbard model. Moreover, since both 
the separation of charge and spin, and the quantum crit- 
ical nature of the spin correlations, are robust features 
of the one dimensional electron gas, they are not limited 
to the integrable Hubbard model, but would be expected 
to survive the inclusion of second and further neighbor 
Coulomb repulsions as well. 
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FIG. 3. The large N = 4n scaling of the pair-binding en- 
ergy, E p (triangles), and spin-gap, A s (circles) of Hubbard 
rings in units of t — 1. All energies are scaled by an explicit 
factor of N to remove the dominant 1/7V dependence, reveal- 
ing the logarithmic dependence predicted in Eq. (^). Solid 
figures are for U — 10 and open figures are for U = 20. 



B. t — J ladders 

Hubbard or t — J chains have no spin-gap in the limit 
A" — > oo, and, in addition, have small superconduct- 
ing susceptibilities irrespective of the doping level. In 
contrast, ladder systems, that is, quasi-one-dimcnsional 
systems of finite width, can exhibit both spin-gap and 
strong tendencies toward superconducting order, even in 
this limit. While these systems are infinite in extent, the 
mesoscopic physics comes in through the finiteness of the 
transverse dimension. Because so much is known theo- 
retically about one-dimensional quantum systems, such 



3 



ladders are also an excellent laboratory for exploring the 
general physical principals set forth in the Introduction. 

The Hamiltonian of a spin- 1/2 Heisenberg ladder, as 
shown in Fig. ||, is 



(0) 



where is a spin 1/2 operator, J is the antiferromag- 
netic exchange interaction, and (ij) signifies nearest- 
neighbor sites, of spacing a, on the ladder. For a Heisen- 
berg ladder, the site occupation numbers, n*, are con- 
stants of motion and are unity. So, the second term 
in the Hamiltonian is simply an irrelevant constant; 
this term will become important for the t — J-ladder 
to be discussed later. For an even-leg ladder of width 
W/a = 2L = 2, the spin-gap is known numerically |23f| 
to be A s w J/2, but this gap vanishes exponentially 
H, as A s - 3.35Jexp[-0.682(W/a)], for wider lad- 
ders, L — 2,3, . . ., as shown in Fig. f|. In contrast, for 
odd-leg ladders, W/a = 2L + 1, Hj has gapless spinon 
excitations, but here, too, there is an important energy 
scale, with the same functional dependence on L as A s , 
below which the physics is dominated by spinons, as for 
a spin 1/2 Heisenberg chain. Since, once again, the spin- 
gap formation is related to superconductivity, as we shall 
show below, this implies that only rather narrow ladders 
(2 or 3-leg) are good candidates for the mesocopic build- 
ing blocks of a high temperature superconductor. 
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FIG. 4. Semi-infinite ladder of coupled quantum S — 1/2 
Heisenberg spins; the width of the ladder is 2L. The plot 
shows that the spin gap A s /J as a function of the width of 
the ladder falls off exponentially fast. 

For a doped system, we can not only ask about the de- 
pendence A. s (x) of the spin-gap and of the pair binding 
energy E p (x) on the concentration of added charge, x, 
but we can also ask about the superconducting suscepti- 
bility. To be concrete, we consider results for the t — J 
model 



Hi 



(T) 



which is defined on the space of no doubly occupied sites, 
and can be viewed as the large U limit of the Hubbard 
model so long as J = At 2 /U <C i. 

Numerical studies ]25j on the two leg t — J ladder with 
J = 0.35t show that the spin-gap is roughly the same 
magnitude as in the undoped two-leg ladder for a range 
of x between and 0.15 and so is the pair-binding en- 
ergy. The three-leg ladder |^6| is more interesting — the 
spin gap is for x < 0.05, then rises to a value compara- 
ble to the spin-gap in the two-leg Heisenberg ladder for 
x ~ 0.15, and then becomes small or possibly zero as x 
gets to be 0.2 or larger. This behavior can be readily 
understood (and was to a large extent anticipated) from 
both perturbative renormalization group analyses [^7|28 
and from strong coupling bosonization |]ll| , ^8| methods. 

Because continuous symmetries cannot be broken in 
one dimension, the ladders can never be truly super- 
conducting. However, their tendency to superconduct 
can clearly be seen in the superconducting susceptibility 
X SC (T); for instance, in a quasi-one dimensional array of 
ladders, one can estimate the superconducting T c from 
the mean-field equation zJxsc{T c ) = 1, where z is the 
number of nearest-neighbor ladders, and J is the inter- 
ladder Josephson coupling. For a one-dimensional sys- 
tem with a spin-gap, the superconducting susceptibility 



Xs 



AT 1/K C 



(8) 



where K c is the charge Luttinger exponent. This expres- 
sion reveals a direct relation between spin-gap and en- 
hanced superconducting fluctuations! A related result, 
which can be more easily compared to the results of nu- 
merical experiments, is the pair-field pair-field correlator, 
D(r), given by 



D{r) ~ A s 



-1/K C 



(9) 



For the purpose of understanding the competition be- 
tween order parameters, discussed below, it is worth not- 
ing that the susceptibility to 2/cf charge-density wave 
(CDW) is 



X.cdw 



A S T 



A',. 



(10) 



<i,j>,cr 



A spin-gap enhances both the CDW and superconducting 
susceptibilities. Indeed, it is tempting, and sometimes 
useful, to think of the spin-gap as a superconducting gap, 
which exists in the absence of superconducting long-range 
order due to violent phase fluctuations characteristic of 
a one dimensional system. But Eq. ( |i"o| ) makes it clear 
that the spin gap can be equally well thought of as a 
CDW gap, where true CDW order has been suppressed 
due to the divergent fluctuations of the acoustic phonons 
in the CDW state. 

The charge Luttinger exponent, K c , is a non- universal 
function of doping, x, and depends on the nature of the 
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interactions. As long as K c > 1/2, the \sc is divergent, 
but not as strongly divergent as the Xcdw, unless K c > 
1. For a Hubbard chain, it is well known that K c < 1, 
but this is not a strict physical bound. For instance, 
by comparing the expression in Eq. (^|) with numerical 
experiments on the 2-leg Hubbard ladder J2{|, one can 
see that K c > 1/2 for a very wide range of parameters, 
and that for some ranges of parameters, K c > 1. For 
instance, for x = 0.0625, K c > 1 for 5 < U/t < 15, that 
is, for intermediate values of U. 



C. Hubbard molecules 

In the past we have examined a number of small Hub- 
bard molecules |32]| by exact diagonalization and found 
pair-binding to operate in all cases. These were the 4x4 
torus, the 8-site cube, and the 12-site truncated tetrahe- 
dron. In all cases, the degeneracy of the state in which 
doping takes place is important, and, in all cases studied, 
second order perturbation theory captures qualitatively 
all the aspects of the phenomenon for what appears to 
be large values of U. This is not so surprising in a finite 
system with finite energy denominators; hence, perturba- 
tion theory is an effective tool in studying large molecules 
such as Ceo, where exact diagonalization methods are 
prohibitive. 

We have also demonstrated that nearest neighbor 
repulsive interactions for the truncated tetrahedron is 
harmful for pair-binding. In contrast, pair-binding in 
t — J ladders remains robust |3(| for nearest neighbor re- 
pulsion, V, as large as 4J. Similarly, realistic frequency 
dependent screened potential for doped Ceo leaves the re- 
sults essentially unchanged from the pure Hubbard model 
|jlf . This we could only demonstrate in second order per- 
turbation theory. 

The conclusion to draw is that in spite of the com- 
plexities at the mesoscopic scale, which surely exist, it 
is well established that certain mesoscopic clusters have 
positive pair-binding energies, E p > 0. Moreover, E p can 
be of order the spin-gap, which is an intrinsically large 
electronic energy scale. There still remains the question 
whether this mesocsopic pairing tendency leads to global 
superconductivity. 



D. The Two-Dimensional t — J model 

Superconductivity in the two dimensional t— J or Hub- 
bard models on a square lattice remains controversial - 
no exact analytic results exist and numerical studies are 
confined to rather small systems or high temperatures. 
However, a few things are clear. The undoped system is 
an ordered antiferromagnet, and therefore A s — 0. More- 
over, while the pair-binding energy for two doped holes is 
positive on finite size molecules, in the physically relevant 



limit of small t/J < 0.4, E p decreases with system size, 
and appears |35| to approach in the thermodynamic 
limit. 

It is plausible, even likely, that on a sufficiently frus- 
trated lattice, the spin 1/2 Heisenberg antiferromagnet 
(the "undoped system" ) can have a quantum disordered 
ground-state with a spin-gap |36|]. Such a system, when 
doped, might also become a high temperature supercon- 
ductor. However, whereas pairing on the mesoscale ap- 
pears to be a very robust phenomenon, finding magnetic 
systems with sufficient frustration to produce a signif- 
icant spin gap is not simple in dimension greater than 
1. So while doping such a homogeneously frustrated an- 
tiferromagnet may provide a mechanism for high tem- 
perature superconductivity, we find pairing on mesoscale 
structures to be a more common (and hence more natu- 
ral) mechanism. 



III. EFFECTIVE INTERACTIONS, COMPETING 
ORDERS, AND SUPERCONDUCTIVITY 

Since a finite cluster cannot be a superconductor, and 
even an extended ladder cannot exhibit a finite temper- 
ature superconducting state, the understanding of this 
mesoscale physics simply serves to define a new, sim- 
pler effective Hamiltonian, which governs the physics on 
longer length scales. 



A. Long-range order in molecular crystals 

When the extended system can be thought of as a 
molecular crystal composed of finite clusters, the effec- 
tive Hamiltonian is obtained by integrating out all the 
molecular orbitals other than those at the Fermi level. 
Where there are several degenerate orbitals at the Fermi 
energy, the resulting effective Hamiltonian will have more 
than one molecular orbital per site. To be concrete, we 
will use as an example the case of a Cgo-like crystal doped 
with x added electrons per molecule, with < x < 6 and 
vary the intra-molecular and the intermolecular couplings 
to examine its properties. 

The effective Hamiltonian is thus of the form 



(11) 



<i]> 



where the effective interaction V on the f h molecule de- 
pends on the total number of added electrons ANj, the 
square of the total spin S| , and the square of the total 
orbital angular momentum, L| (or more precisely, the 
equivalent Casimirs of the ichosohedral group), and Hij 
is the intermolecular interaction. In the simplest case, 
Hij can be taken to consist of a sum of intermolecu- 
lar hopping operators, which transfer a single electron 
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between molecules i and j, but it can also contain ad- 
ditional terms, such as assisted hopping and Coulomb 
interactions. 

Where E p is positive on the cluster, V can be thought 
of as an effective attraction between pairs of electrons 
on a given molecule. The ratio of E p to the intermolec- 
ular band- width, W, is an important parameter in the 
problem. 

1. Strong- coupling limit 

If the typical magnitude of V is large compared to 
W, we can further thin the low-energy Hilbert space 
by integrating out all states other than the molecular 
ground-states in each charge sector. In most cases, this 
limit is unphysical, but because it is easily analyzed, it is 
worth considering for pedagogic purposes. In particular, 
if E p > 0, we can eliminate all states with an odd number 
of electrons per site, while for E p < 0, we can eliminate 
all states with an even number of electrons. The resulting 
effective Hamiltonian has no remaining Fermionic degrees 
of freedom; the states are specified entirely by the num- 
ber of electrons, and possibly the spin, and the orbital 
moment of each molecule. Under appropriate circum- 
stances, we are led to spin or orbital effective Heisenberg 
models, and even to a bosonic version of the t — J model. 

Depending on details of the interactions, such effec- 
tive Hamiltonians can easily lead to a rich variety of 
ordered phases, including fcrro or antiferromagnets, or- 
bital ferromagnets or antiferromagnets, spin or orbital 
nematics, modulated (incommensurate or high-order 
commensurate) versions of any of these states coexisting 
with charge-density wave order, etc. This is, of course, 
both exciting and discouraging. It suggests that in the 
presence of strong interactions in systems with mesoscale 
structure, there should appear an enormous variety of or- 
dered phases, some of them are new states of matter that 
have never before been documented. While it invites us 
to exercise our imagination, the prospects of predicting 
the ultimate form of order from microscopic considera- 
tions in any such system becomes a difficult task, if not 
impossible. 

To take the simplest case, let us imagine that the ef- 
fective interactions favor an even number of electrons per 
molecule (Ep > 0), and disfavor molecular states with 
non-zero spin or orbital moment. In this case, the low 
energy physics of the molecular crystal reduces to that 
of interacting bosons (electron pairs) on the molecular 
lattice 

ir Jf 5>*&& • v(pfa)] 

3 

£ /,,>:, '< a • '>:'>, • .... (12) 

<ij> 



Here, is a charge 2e bosonic creation operator, /x& is 
1/2 the electron chemical potential, V is the repulsion 
between pairs on the same molecule and includes a hard- 
core interaction that forbids more than three bosons on 
the same molecule, t pa i r is the pair hopping term which is 
second order in the electron hopping matrix element, and 
. . . signifies additional interactions between bosons on 
neighboring molecules. Of course, even here, the ground- 
state phases of the system will depend on the precise bal- 
ance between these residual interactions, as well as on the 
mean boson density. 

Typically, if the electronic charge per molecule is x = 2 
or 4, the ground-state will be insulating, with essentially 
a fixed number of electron pairs on each molecule. This 
state has no broken symmetry, but is a Mott insulator, 
in the sense that in the non-interacting limit the system 
would be metallic. For other mean-charge densities, this 
system will either form a pair-density wave (presumably 
insulating) or a superfluid. In the latter case, the pair- 
ing scale is intramolecular, and hence large, while the 
superfluid density, and hence the global phase ordering 
temperature, T c , will be smaller, determined by the in- 
termolecular bandwidth. Moreover, in this case, pairing 
should produce a pseudo-gap which survives at temper- 
atures well above T c . In this limit, because the effective 
boson kinetic energy is proportional to the square of the 
inter-molecular hopping matrix elements, one would ex- 
pect T c to be a decreasing function of lattice parameter, 
or an increasing function of pressure. 

2. Weak coupling limit 

If the typical scale of interactions, V, is small compared 
to W, and if the molecular crystal is three dimensional, 
it is reasonable to treat the effective Hamiltonian in a 
mean- field approximation, that is, by variational ground- 
states, in the spirit of the BCS state, by solving the non- 
interacting problem in the presence of a symmetry break- 
ing mean-field. Because of the narrowness of the bands 
in typical molecular solids, it is equally unclear whether 
this limit is actually realized in real materials. But it is 
at least another limit in which controlled theoretical re- 
sults can be obtained. It is probably closer to reality than 
the strong coupling limit, at least in relatively metallic 
materials, such as alkali doped Cgo- 

In this limit, the various charge, spin, and orbital mo- 
ment density wave states are unlikely, unless there are 
particular nesting wave- vectors of the relevant Fermi sur- 
face. Thus, except at certain commensurate values of x, 
they can be safely ignored. 

For positive E p , a BCS instability to a superconducting 
state is still generic with a superconducting transition 
temperature, 

T c ~ We(V cS ) exp[-l/p(E F )V cS (x)} (13) 
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where is the Heaviside function, and 

p{E F )(xl/W (14) 

is the density of states at the Fermi energy of the non- 
interacting band-structure. The quantity 

V cS {x) = p x(Q, S)E P (Q, S) (15) 
Q,S 

is the average of the pair-binding energy for the various 
charge states of the molecule averaged over the appro- 
priate probability distribution (J2q s P x (Q, S) — 1) of 
charge and spin-states of the molecule. Here, T c would 
be expected to be a strongly increasing function of lat- 
tice parameter (decreasing function of pressure) since it 
varies with the inverse exponential of l/W. 

While it might seem at first that this BCS expression 
for T c is generic to any mechanism of superconductiv- 
ity with an intra-molecular pairing force (for instance, 
one generated by coupling to the intra-molecular optical 
phonons), the strong x dependence of V implied by the 
electronic mechanism discussed here is very pronounced 
and unique. For instance, from the perturbative results 
for the Hubbard model on a Ceo cluster described above 
||, one might expect that for intermediate £/, V cS is pos- 
itive (attractive) for x — 3, and negative (repulsive) for 
x — 2 or x = 4. Crudely, this leads to a doping dependent 
effective attraction of the general form 

V cS = V (x-x )(x 1 -x), (16) 

with 2 < xo < 3 < x\ < 4. This leads to a large T c only 
in a constrained region about x = 3, and to T c = for 
x < xq or x > x\. This is illustrated in Fig. 5, where, 
to make a crude comparison with experiments, we have 
taken W = O.leV, p(E F )V = 1/3, x Q = 2, and x t = 4. 

T C (K) 




2.5 3 3.5 4 

X 

FIG. 5. T c as a function of electron concentration x, for the 
simple, weak-coupling model of correlation-induced pairing in 
Ceo described in Eqs. @ and (Q. 



B. Long-range order in arrays of ladders 

Long-range order in an array of ladders is related to the 
general problem of order in quasi-one dimensional elec- 
tronic systems and has been studied for for many years 
; this is theoretically one of the best understood cases 
of competing orders. In particular, the occurrence of a 
spin gap enhances the tendency for both CDW and su- 
perconducting order in a way that does not distinguish 
between them. 

Even in this case, the competition plays out in a man- 
ner that was not previously recognized. It is not suffi- 
cient to simply identify the most divergent susceptibility 
of the isolated ladder in order to determine which type 
of long-range order will win out. There are, in addi- 
tion, a variety of marginal interladder forward scatter- 
ing interactions which can affect the balance, or whether 
or not there is an ordered ground-state at all |3. A 
cautious, but optimistic conclusion from the mesoscale 
calculations is that the spin-gap, which can be large on 
narrow ladders, permits the possibility of high tempera- 
ture superconductivity, but it can also lead to a variety 
of alternative insulating or metallic phases. 

IV. FINAL REMARKS 

Superconductivity in simple metals is well understood 
in terms of the BCS mechanism — the normal state is 
well approximated as a Fermi liquid, and the effective at- 
traction between electrons is mediated by phonons. The 
effects of electron-electron interactions appear in this de- 
scription in the form of a few Fermi liquid parameters 
(which are typically small), and in the Coulomb pseudo- 
potential, /!*. The Fermi liquid character of simple met- 
als is exceedingly robust - in the absence of special nesting 
conditions of the Fermi surface or very strong interac- 
tions, the only instability of the metallic state is the BCS 
instability. A conventional superconductor, when driven 
normal by suitable alloying, by weak disorder (magnetic 
impurities), increasing temperature, or magnetic field, re- 
mains a highly conducting metal. 

In both the cuprates and the organic superconductors, 
including doped Cgo, there is strong reason to doubt 
the accuracy of the Fermi liquid description of the nor- 
mal state. Screening is generally poor, and the typical 
magnitude of Coulomb energies is large compared to the 
valence-electron bandwidth. And, indeed, the materi- 
als are not robustly metallic - resistances greater than 
the quantum of resistance (often accompanied by "insu- 
lating" negative values of dp/dT) are frequently found 
in these materials when the superconducting state is de- 
stroyed. Even when the resistance has a metallic temper- 
ature dependence (positive dp/dT), as it does in doped 
Ceo in the "normal" state, the resistance is too large 
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to correspond to any sort of freely propagating quasi- 
particle. In electron-doped Cgo, the room temperature 
resistance |58]jll| is 2-5 mfi-cm, which were it interpreted 
in terms of a Drude conductivity would correspond to a 
mean- free path of 1-2 A; in hole doped Cgo, the room 
temperature resistance |l3j is almost 20mf2-cm, which at 
face value would correspond to a mean-free path of or- 
der O.lAl Finally, in addition to superconducting order, 
evidence of a variety of competing orders abounds. 

Given the striking evidence of strong correlation ef- 
fects in these materials, we suppose that they are funda- 
mental to the mechanism of their unprecedentedly high 
T c 's. However, we suggest that such a mechanism is ro- 
bust only if the pairing originates on mesoscale struc- 
tures. Therefore, it is not accidental that the materials 
in question are either molecular crystals, with reason- 
ably large molecular building-blocks, or materials which 
exhibit self-organized mesoscale structures. 



A. The electron-phonon coupling 

Electron-phonon interactions are reasonably strong in 
all the high temperature superconductors. In particular, 
the phonon frequencies in Cqq are much higher than in 
conventional metals. In spite of this, we very much doubt 
that superconductivity is driven principally by electron- 
phonon interactions. That T c exhibits isotope effect can- 
not be an argument as it can be explained by the elec- 
tronic mechanism as well p£| . 

In the first place, the pairing induced by strong repul- 
sive interactions on all the clusters that have been stud- 
ied to date have d-wave character, for reasons that are by 
now well understood. Phonons typically produce s-wave 
attraction, and are often pair-breaking in other channels. 
Indeed, when we jio) analyzed the effects of electron- 
phonon couplings on an assumed electronic mechanism of 
pairing on a Ceo molecule, we found that most, but not 
all of the high energy phonons are pair-breaking. Note 
that local d-wave character may or may not imply that 
the global order has d-wave symmetry — for instance, 
we have found [mj^] that the superconducting order in 
electron-doped Cgo is globally s-wave. 

In the second place, the strong electron-phonon inter- 
action, presumably responsible for high T c , is always ac- 
companied by self-trapping, or bipolaron formation. The 
result is a large exponential Frank-Condon reduction of 
the energy scale for coherent motion of charge, inevitably 
leading to an insulating state. 

In the weak-coupling limit, the usual BCS result ap- 
plies with an exponentially small pairing scale, and the 
retarded nature of the electron-phonon coupling plays an 
essential role — the Coulomb interaction between elec- 
trons is reduced due to the fact that they do not have 
to be at the same place at the same time. However, in 
the cuprates, and even more clearly in doped Cgo, the 



phonon frequencies are comparable to the electron band- 
width, so retardation does not lead [|d| to any reduction 

of (J,*. 



B. Superconductivity in the polyacenes 

Superconductivity has been recently discovered in a 
surface layer of electrostatically doped crystals of an- 
thracene, tetracene, and pentacene [jl4|. These molecules 
can be viewed as segments of a two-leg ladder with every 
other rung weakened. Thus, we believe that the basic 
mechanism of superconducting pairing is the same as the 
one discussed previously for the two leg t — J and Hub- 
bard ladders. As a consequence, we are lead to expect 
that the superconducting order will have a d-wave-like 
symmetry, the order parameter changing sign under 90° 
rotation. However, because of the quasi-one dimensional 
geometry, it is possible that there are no gapless "nodal" 
quasi-particles in this system, despite this [42). 

We note that a number of features of the experiment 
are consistent with this suggestion. In the first place, T c 
is maximal when there is approximately an odd number 
of electrons per molecule. Moreover, the fact that the op- 
timal T c decreases with the increasing size of the molecule 
is consistent with a mesoscopic origin of the pairing. 



C. C60 

Narrow bands in Cqq make superconductivity at 
high temperatures (presently 52 K) difficult unless the 
molecule itself is in some sense superconducting. The 
same difficulty holds for all narrow band materials f43j ] 
for the following reasons. 

First, the net interaction at u> = must be repulsive, 
otherwise we will get a CDW, that is, the interaction is 



V cS (q,uj) 



Aire 1 



q 2 e(q,u) 



e{q,LU = 0) > 0. 



(17) 



Typically, a peaked density of states implies that only 
a single band is relevant, corresponding precisely to the 
band in which the peak in the density of states is situated. 
Thus, if only one band is relevant, no local held effects are 
possible, and the ions will see the same e as the electrons. 

The second reason is that the retardation of the 
electron-phonon interaction in BCS superconductors 
(that is so crucial for the existence of a net electron- 
electron attraction through the reduction of the Coulomb 
pscudopotential) is no longer operative. This is because 
the Coulomb pseudopoetntial (i* is given by 



1 + M ln^ ' 



(18) 



where fj, is the Coulomb matrix element at the bare high 
energy scale, W is the band width, and luq is the typical 
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phonon scale. But if W is not much larger than u>o, which 
is manifestly the case in Ceo, we get little reduction in 
H*. 

In contrast, in a molecular crystal, the existence of a 
structure at an intermediate scale implies atomic energy 
scales much larger than intramolecular energy scales, 
which in turn are larger than the intermolecular scales 
of the solid. This clear separation of scales allows us 
to solve first the problem of a single Cgo molecule alone 
and then incorporate the coupling between the molecules. 
This approach leads to the by now well substantiated 
phenomenon of attraction from repulsion that has been 
the corner stone of our approach. Even issues concern- 
ing the proper screening of the longer range pieces of the 
Coulomb interaction can be reliably analyzed |n| taking 
advantage of this hierarchy of energy scales. 

We conclude with a few observations regarding the 
present experiments fl2] , |l3|| on Geo FET's. It is clear 
that high density of states cannot be the principal fea- 
ture responsible for the high temperature superconduc- 
tivity in hole doped Ceo as T c is clearly not peaked at 
the same charge per molecule as the density of states. 
Moreover, the T c profile as a function of doping does not 
presently have any natural explanation in the electron- 
phonon theory of a BCS superconductor; there appears 
to be no reason for the electron-phonon coupling to show 
such a strong dependence on the band filling. In the 
electron doped case |E3], the onset of superconductivity 
is very sharp. This is consistent with our theory in which 
T c is highest at a doping corresponding to an odd number 
of added electrons or holes, where the pair-binding on a 
single molecule is largest, but then vanishes for doping 
corresponding to an even number of added electrons or 
holes [||. A crude model was shown in Fig.||. Of course, 
sample inhomogeneity and disorder can broaden the dop- 
ing profile, which one should be able to better address in 
the near future as the sample quality improves further. 
Moreover, relatively low temperature tails in T c vs x could 
result from residual electron-phonon interactions. 

If the explanation of superconductivity lies in the elec- 
tronic correlation effects on the mesoscopic scale, there is 
no reason why T c could not be raised further (> 52 K) in 
the near future. Such high values of T c will constitute a 
strong argument against the electron-phonon mechanism 
of superconductivity in these materials. 
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APPENDIX A: PERTURBATION THEORY FOR 
HUBBARD CHAIN 

First order perturbation theory leads to the conclusion, 
which turns out to be exact for all U/t, that for 
N = An, the ground-state is a spin singlet for Q = and 
Q = ±2 and a doublet for Q = ±1, while for N = An + 2 
the ground-state is a singlet for Q — 0, a doublet for 
Q = ±1, and a triplet for Q — ±2. As a consequence, 

E P {V) = ^^ + 0{UVt*) (Al) 

where Ai, a number of order unity, is positive for N — An 
and negative for N = An + 2, reflecting the difference in 
spin of the doubly charged state. 

For N = An + 2 and U = 0, the ground-state with 
Q = is a unique spin singlet with a gap; for Q = ±1, 
the ground-state is still unique up to symmetry — it has 
spin 1/2 and non-zero crystal momentum; however, for 
Q = ±2, there are three singlet and one spin triplet states 
that are degenerate. In first order in U/t, the degeneracy 
of the doubly charged system is lifted and, as required by 
Hund's rule, the spin-1 state is selected as the ground- 
state. In the Hubbard model, there is no direct interac- 
tion between electrons of like spin, so a consequence of 
the triplet character of the doubly charged state is that 
there is no first order contribution to the pair-binding 
energy. However, mainly because of the unfavorable fer- 
romagnetic correlations in the doubly charged state, it 
turns out that A^{An + 2) is negative — there is pair 
repulsion. 

For N — An, there are two degenerate orbitals at the 
Fermi level of the non-interacting neutral molecule re- 
sulting in accidental degeneracies. To first order in U/t, 
the energies of two of the singlet states of the neutral 
molecule are elevated, but one singlet state remains de- 
generate with the triplet state. This miracle is due to 
an Umklapp process which allows a singlet pair of elec- 
trons, each with crystal momentum 7r/2, to scatter into 
the state in which each electron has crystal momentum 
— 7r/2. The final degeneracy for Q = is then lifted in 
second order perturbation theory, in which the antifer- 
romagnetic correlations of the singlet state are preferred 
over the ferromagnetic (triplet) state. This violation of 
Hund's rule is peculiar to the case of N — An and Q = 0. 
By the same token, the antiferromagnetic correlations of 
the neutral molecule lead to anomalous stability of the 
neutral state, and hence a negative A2(An) — there is 
pair binding. It is worth noting that, although for large 
U/t the perturbative results are, of course, quantitatively 
poor, as there are no groundstate level crossings as a 
function of U/t, perturbation theory gives correct results 
for the ground-state quantum numbers of all these states. 
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